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Abstract. The effect of sequence heterogeneity on the dynamics of molecular motors 
is reviewed and analyzed using a set of recently introduced lattice models. First, we 
review results for the influence of heterogenous tracks such as a single-strand of DNA 
or RNA on the dynamics of the motors. We stress how the predicted behavior might be 
observed experimentally in anomalous drift and diffusion of motors over a wide range of 
parameters near the stall force and discuss the extreme limit of strongly biased motors 
with one-way hopping. We then consider the dynamics in an environment containing 
a variety of different fuels which supply chemical energy for the motor motion, cither 
on a heterogeneous or on a periodic track. The results for motion along a periodic 
track arc relevant to kincsin motors in a solution with a mixture of different nucleotide 
triphosphate fuel sources. 
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1. Introduction 

The study of molecular motors has been transformed in recent years with the increasing 
use of single molecule experiments P3 |2j. In one key experiment an external force is 
applied to a molecular motor opposing its motion [21 El El Hj • Typically, as the force is 
increased, the velocity of the motor decreases until it is completely stalled. The behavior 
of the velocity as a function of force provides much information on the chemical cycle 
underlying the motion of the motor. For example, the stall force is a direct estimate of 
the force exerted by the molecular motor. The experiments have also motivated much 
theoretical work on the dynamics of the motors [3 IHl El EI] • Fitting the experimentally 
obtained velocity-force curves allows extraction of detailed information on the chemical 
cycle of the motor [TT1IT21IT3]. 

Most theoretical studies have focused on motors which move on featureless, or 
periodic linear tracks [3 El El HH EE3 EH- Such a description would be appropriate, 
for example, for kinesin which moves along a microtubule filament, which is periodic, 
using only ATP for its motion However, in many cases the assumption of a 

periodic medium fails. Examples of motors which move on heterogeneous tracks include 
RNA polymerase [U El which moves along DNA, ribosomes which move along mRNA, 
helicases [H3 El which unwind DNA, exonucleases [T51 IT§] which turn double-stranded 
DNA into single-stranded DNA and many others. All these motors move along tracks 
which are inherently "disordered" or heterogeneous due to the underlying sequence of the 
linear template. Theoretically, molecular motors moving along disordered tracks have 
received much less attention [201 [^1 [22]- Another form of heterogeneity which has largely 
been ignored arises from the different chemical fuels which may be used by molecular 
motors to move along the track. For example, RNA polymerase uses different nucleotide 
triphosphates (NTP's) which build the mRNA it produces, each supplying a different 
amount of chemical energy, to move along a DNA strand. A different "annealed" form 
of disorder (in contrast to the "quenched" disorder embodied in a particular nucleotide 
sequence) can be present even in molecular motors moving along perfectly periodic 
tracks in a solution containing several distinct types of chemical fuels. For example, 
it is known that kinesin can move using other nucleotide triphosphates (such as GTP) 
instead of ATP, albeit less efficiently [231 EH EH]- 

Recently, we have introduced a simplified model for molecular motors which allows 
the effects of disorder to be studied in considerable detail [25 • We have focused so far on 
heterogeneous tracks and argued that near the stall force the dynamics of the motor is 
strongly affected by the heterogeneity embodied in a particular DNA or RNA sequence. 
Due to the "sequence disorder" on which the motor is moving (many DNA sequences 
have only short range correlations j2E]) the displacement of the motor as a function 
of time ceases to be linear in time close enough to the stall force. The displacement 
becomes sublinear in time, growing as t M , with fi varying continuously from 1 to as the 
stall force is approached. As discussed below there are also anomalies in the diffusive 
spreading about the average motor position which extend even further below the stall 
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In this paper we review some of these results, stressing several experiments which 
could be performed to test the predictions of the model. We also explore the effect 
of heterogeneous fuels on the motion of molecular motors. In [2T] it was suggested 
that inhomogeneous fuel concentrations could enhance significantly the regime near the 
stall force over which anomalous dynamics is observed. Here we study this type of 
disorder numerically and illustrate the dramatic effect of varying the concentrations of 
the different fuels used to power the motor. For motors moving along heterogeneous 
tracks we also discuss the dynamics in the extreme limit where detailed balance is 
violated and motors never take backward steps. Finally, we consider motors moving 
along a periodic substrate powered by different kinds of fuels. We discuss, for simple 
cases, the expected velocity of the motor as the relative proportion of two different 
types of fuel in the solution is varied. The behavior of more complicated models is also 
discussed. 

2. The Model 

In this section we define the model used throughout the paper. We start with a 
special case of the general class of n-state models explored by Kolomeisky and Fisher 
[TT| IT2"] IT3*1 l2"Tj and consider a "minimal" motor with only two internal states. This 
simplified model reproduces important features of previously studied systems and allows 
us to explore generic behavior in new situations in a minimal form. The model is easily 
generalized to account for heterogeneous fuels and tracks. When appropriate, we will 
mention how results are modified for general n-state models. The model has been 
introduced and studied in detail in |2^ and here we only review its basic properties. 
We begin by assuming a perfectly periodic substrate. The location along the one- 
dimensional track, x, is assumed to take a discrete set of values x m , where m — 0, 1, 2 . . . 
labels distinct a and b sites. Although not essential, we assume for simplicity that the 
distances between x m+ i — x m and x m+2 — £ m +i are equal and set x m+2 — x m = 2a , 
which is the size of a step taken by the motor after completing a chemical cycle such as 
hydrolysis of ATP. In general, as discussed in the distance traveled by the motor 
between internal states may be different for different internal transitions. However, we do 
not expect such modifications to affect the long time behavior over a range of parameters 
near the stall force. The dynamics embodied in the model is shown schematically in 
Fig. d Internal states labeled by a have an energy e = while internal states labeled 
by b have a higher energy e = As. The local detailed balance condition (in temperature 
units such that ks = 1) is satisfied by our choice of rate constants, 
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Figure 1. Graphical representation of a simplified lattice model for molecular motors 



and the relevant rates wl* , 



and energy difference Ae. The distinct even 



and odd sublattices are denoted by a and b respectively. 



Following Ref. [7], there are two parallel channels for the motion. The first, represented 
by contributions containing a and a', arise from utilization of chemical energy biased 
by a chemical potential difference A/x between, say ATP and the products of hydrolysis 
ADP and Pj. The second channel, represented by the terms containing uj and uj', 
correspond to thermal transitions unassisted by the chemical energy. We assume that 
the externally applied force F biases the motion in a particularly simple way (consistent 
with detailed balance) and define / = Fa®. If the substrate lacks inversion symmetry 
(a necessary condition for directed motion, driven by A/x, when / = [7j), we have 
a' a and uj' ^ uj. If the fuel is ATP, the chemical potential difference which drives 
the motion is ^1] 

*- r Ki&)-Kp§ik)]- (2) 

where the square brackets [. . .] denote concentrations under experimental conditions and 



the brackets 



J cq 



denote the corresponding concentrations at equilibrium. 



The rate constants in Eq. Q define a set of differential equations for the probability 
P n (t) of being at site n at time t. For odd n one has 
dP n (t) 



dt 

while for even n 
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It is illuminating, especially when we consider rate constants which depend on the 
position along a heterogenous track, to study two limits of these equations. In the 
first the chemical potential difference Afi and the applied force / are small compared 
to the energy difference Ae so that b states relax quickly compared to a states. This 
condition implies that (w^* + w£~) ^> (w^ + w^) so that in the long time-limit to a 
good approximation the left hand side of Eq. El may be set to zero. Upon solving for 
P n (t) with n odd and substituting into Eq. HI we obtain differential equations just for 
the even sites 
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(n even) . (5) 

Similarly in the limit All Ae (with / near the stall force) the motor spends most if 
its time in b states. Now, in the long time-limit to a good approximation the left hand 
side of Eq. 0] may be set to zero. The remaining differential equations for the odd sites 
read 

dP n (t) = w^w^P n+2 (t) + w^w^P n - 2 (t) - (w?w? + wfw^)P n (t) 
dt (w^ + wf) 

(n odd) . (6) 

Note that in both limits the dynamics of the motors on long-times can be described 
by a random walker moving on an effective energy landscape associated with what 
is in general a non-equilibrium dynamics. Upon absorbing the denominator factors 
(w^ + w'a) and {w^ + into a rescaling of the rate constants in the numerator, 
the effective energy landscape can be read off from Eq. El or Eq. El One finds that the 
effective energy difference between two sites which are two monomers apart is given by 

E n+2 — E n = AE = T In (^^) . (7) 

For a periodic track, this leads to a tilted energy landscape (with tilt controlled by 
A [a and /) and an effective energy difference between 2m-adjacent monomers ImAE. 
The tilted energy landscape leads to diffusion with drift on long time-scales and large 
length-scales. The effective energy landscape can also be obtained by assuming detailed 
balance and equating the rate asymmetry between two neighboring even sites to an 
effective energy difference AE. It is straightforward to verify from Eq. [7] with Eq. ^ 
that for a periodic substrate no net motion is generated when the external force / = 
and the chemical potential difference All = 0. Also, when there is directional symmetry 
in the transition rates, a = a', u — lu', and / = no net motion is generated even when 
All ^ 0. Absent this symmetry, chemical energy can be converted to motion. These 
conditions are equivalent to those presented in EI| for continuum models and are 
exhibited here in a minimal model. The effect of the externally applied force is simply 
to bias the motion of the motor in the direction in which it is applied. 

The velocity for a motor moving along a periodic track in the two limits discussed 
above can be obtained by taking the continuum limit and yields v = (w^w^* — 
w^w^~) I (w^ + w^~) for the limit specified by Eq. Eland v = {w^w^ — w^w^~) / \w~^ + 
w^~) for the limit specified by Eq. El More generally, for periodic rates, it is 
straightforward to calcualte the velocity, for example using Bloch eigenf unctions \k) ~ 
e %kx anc j expanding the eigenvalues in the wavevector k. The linear term gives the 
velocity and the quadratic part the effective diffusion constant. For the velocity one 
finds 

v = - (8) 

^ + *v + w t + w r 

An alternative to studying the solution of the equations, which will be very useful 
throughout this paper, is to use Monte-Carlo simulations. The rates specified in Eq. 
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(JTJ can be simulated using the following procedure: To make the simulation efficient we 
first normalize the entering or leaving rates for a site so that the largest one is unity. 
Then, at each step we choose with equal probability attempting to move the motor 
to the right or left on the lattice. Following this choice a random number is drawn 
from a uniform distribution in the interval [0,1]. The motor is moved in the chosen 
direction provided the random number is smaller than the corresponding rate. Thus, if 
a motor finds itself on site a in Fig. 1 and > w^, w^, is the largest rate it will 
(after the rescaling) move one step to the right with probability 1/2, one step to the left 
with probability 1/2 (w^/w^*) and it will stay put with probability 1/2[1 — (w^~ /w^)]. 
Note that the probability of actually moving one step (right or left) to the 6-sublattice 
during a particular attempt is 1/2[1 + {w^ /w^)}. Once the 6-sublattice is reached the 
procedure is repeated with the rates and (note that the probabilities are still 
obtained by dividing by w^). To compare, for example, velocities for different choices 
of rates the overall number of attempts is rescaled at the end by the fastest rate (taken 
to be in the above example). The same procedure is followed for both homogeneous 
and heterogeneous tracks. For the latter the largest rate is chosen from all possible 
hopping rates along the track. This protocol ensures relaxation to equilibrium in the 
absence of chemical or mechanics driving forces |2*H] . 

In [21] we have analyzed in detail the motion of the model when the track is not 
periodic. Such tracks arise naturally, for example, for motors such as RNA polymerase 
or helicases which move on DNA which has a well defined sequence. In this case the 
energy difference AE(m) now becomes an explicit function of the location m along the 
track, due to the dependance of the rates on the location on the track. 

To understand the energy which arises for heterogeneous tracks consider the 
"integrating out" procedure applied to the three sites shown in Fig. |2l where three 
distinct motor binding energies, Ei, E 2 and E 3 , are indicated explicitly. We work in the 
limit A/i 3> AE, with AE = E\ — E% or AE = E3 — E%, and / close to the stall force, 
so that the approximation leading to Eq. El ( "integrating out" site 2) is appropriate. As 
rates for the heterogeneous cluster shown in Fig. we take 



where the arguments "(13)" appended to the u>'s, a's,a"s, u/s and a/'s simply mean 
that these are the heterogenous rates appropriate to the cluster 1 — 2 — 3. These rates 
obey detailed balance conditions for the two channels, and have a similar dependence 
on Ayu(13) and / and various energy differences as the rates in Eq. ^ The notation 
A/i(13) indicates that the chemical potential difference could depend on which NTP 
(in the case of RNA polymerase) provides the energy for that particular step. Upon 



«C(13) 
<(13) 

«C(13) 
<(13) 



[ a (13) e M13)/T + w(13)]e -( B2 - El )/T-//2T 

[a(13) + cj(13)]e //2T 

[ a '(13) e A*13)/T + ul{13)]e - { E 2 -E 3 )/T + f/2T 

[a'(13)+u'(13)]e- f/2T . 



(9) 
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Figure 2. The transition rates and energy levels of three sites, spanning two monomers 
corresponding to a motor moving along a heterogenous track. Note that E\ ^ £"3 due 
to the non-periodicity of the track. The transition rates now depend explicitly on the 
location along the track. To avoid cluttering this dependence was suppressed in the 
figure. 



assuming fast relaxation of site 2 in Fig. |2l from a formula similar to Eq. [7[ 

AF -F F + 2f Tl J Kl3)e^ + ^(13))K(13) + ^(13)) - 
AE 13 -E 3 -E 1 + 2f-T\n [ (Q , (13)e A^ + w / (13))(a(13) + w(13)) _ 

= E 3 -E 1 + 2f + r) 13 (10) 

Eq. El illustrates the following important points, applicable to motor molecules on 
heterogenous tracks more generally: (a) if AyU,(13) = 0, then E\ 3 = E 3 — Ei + 2f, with a 
similar formula for all neighboring pairs of odd sites. Thus, in the absence of chemical 
energy, we have a "random energy landscape" with bounded energy fluctuations; (b) 
if a(13) = a'(13) and u;(13) = u/(13) (inversion symmetry), chemical energy does 
not lead to net motion between sites 1 and 3, as discussed above; (c) in general, 
AE V3 = E 3 — Ex + 2f + rj 13 , where 7713 is a random function of position along the 
heterogeneous track. Upon passing to a coarse-grained position 77(771), where m is the 
position along the track we see that the effective 'coarse grained' energy difference 
between two points m\ and m-i = mi + m which are m-monomers apart is given by 
77(772'). The landscape itself behaves like a random walk with fluctuations which 
grow as v^m, corresponding to a random-forcing energy landscape (for RNA polymerase, 
the different chemical potentials of the nucleotides in the transcript also contribute to 
a random forcing landscape |21j). 

The self-similar structure of the random force landscape leads to interesting 
dynamics near the stall force. As the stall force is approached the dynamics slows 
down and becomes dominated by motion between deep minima of the energy landscape 
|29j . The minima correspond to specific locations along the track where the motor 
tends to pause. The distribution of dwell times at these minima, P(t), averaged over 
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the different locations on the track, is expected to behave as r~( 1+At ), where fi (not 
to be confused with a chemical potential!) is related to the force, fluctuations in the 
effective energy landscape and temperature. For random forcing energy landscape where 
the energy difference between two points is drawn from a Gaussian distribution with a 
variance V = r](m) 2 , where the overline denotes an average along the sequence, one can 
show that HH 

//(/) = 2T|A£ /=0 - 2f\/V , (11) 

where AEj = q is the mean slope of the potential (averaged along the sequence) at zero 
force. The exponent fi thus decreases continuously to zero as / increases toward the 
stall force of the motor (defined by /x(/ s ) = 0). For more general distributions of the 
effective energy difference the value of \i might be different from Eq. ^2 by factors of 
order unity. 

Near the stall force the expected distribution of pause times becomes broader as \x 
becomes closer to 0. The dynamics of the motors are altered from diffusion with drift 
when the pause-time distribution becomes very broad. The dynamics then depends on 
the numerical value of \i defined in Eq. ^2 [2H 1221 EH] : 

• \x < 1 - Around the stall force, both the drift and diffusive behavior of the motor 
become anomalous. The displacement of the motor as a function of time increases 
as t M . Thus, in this region the velocity is undefined, in the sense that it depends 
on the experimental observation time, through v ~ t^T 1 . Moreover, the spread 
of the probability distribution of the motor about its mean position also behaves 
anomalously with a variance which grows as t£. Experimentally, for a given 
this anomaly should lead to a convex velocity as a function of force curve in the 
vicinity of the stall force. The curve will become more and more convex as t^ is 
increased; the velocity actually vanishes for a range of /'s near the stall force in 
the limit — > oo. 

• 1 < fi < 2 - Further away from the stall force the displacement of the motor as a 
function of time grows linearly. At long-times the velocity becomes independent of 
the averaging window. However, the variance of the probability distribution around 
the mean is anomalous and grows as t 2 ^ . 

• fj, > 2 - Far below the stall force both the displacement and the variance of the 
probability distribution around the mean grow linearly in time, as in conventional 
diffusion with drift. 

These results can easily be shown to apply as well to general n-state models. 
Moreover, it can be argued that even if several parallel channels exist for moving from 
one monomer to another the results are also qualitatively unchanged [21]. 

Experimentally, the predictions of the model can be tested by measuring the 
displacement of the motor as a function of time, averaged over different experimental 
runs (and, possibly, sequences). Each time trace of the motor position will have an 
irregular shape due to pauses, which will increase in duration as the stall force is 
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approached, at specific locations along the track. However, averaged over many time 
traces (or sequences) the expected displacement will grows as V 1 with \i < 1 close enough 
to the stall force. Note that if one averages over time traces for a fixed sequence, 
the displacement is expected to grow as where s(t) has fluctuations of order 

unity, because the sequence information is not completely erased in this case. An 
alternative experimental test would be to measure the distribution of dwell times V[f). 
Because V{j) ~ l/r 1+fl for large r, the distribution becomes wider as the stall force is 
approached. Monitoring V(t) has the advantage of probing the wide distribution even 
in regimes which are not very close to the stall force. 

We now consider limiting cases of the above model on heterogeneous tracks. These 
illustrate a number of interesting features and suggest ways in which the anomalous 
dynamics might be observed experimentally. We will also use the model to explore 
heterogeneous chemical energy sources for motors on a periodic track. This situation 
may be realized in motors such as kinesin which can use several types of chemical 
energy to move along the track. From the two state model described above we deduce 
the expected behavior of the velocity as the relative proportions of the different fuels is 
varied and mention generalizations to more general n-state models. 

3. Strongly biased motors 

In this section we study motors moving on a heterogeneous substrate in the limit where 
one of the transition rates is strongly biased in a certain direction. An extreme limit 
occurs when one of the transition rates, in, say, the backward direction, is zero. Although 
this limit violates detailed balance, it could be a reasonable approximation for certain 
strongly biased experiments. One such model is a special case of the two state model 
discussed in Section El 

w ? = {ae^' T + uj)e-^' T ~fl 2T 
w f = ( a + uj)e f/2T 

«£- = J e -MT + f/2T (12) 

-> / -f/2T 

w b = uo e J 1 

Here we have set a' = so that in the A/j, 3> T limit the motor will be strongly biased 
to move towards the right. Physically, this situation corresponds to a motor which can 
use chemical energy only to move in a certain direction. (We expect qualitatively similar 
results for a wide variety of strongly biased "one way" models). Next, we assume an 
extremely strong bias A/x 2> T limit of the model such that ae A/X//T is very large but 
a and u are so small that can be set to be zero. This limit only makes sense far 
from the stall force. The stall force of a model with wl~ = (as any other model were 
one of the reactions is assumed to be unidirectional) is infinite: The effective energy 
landscape, Eq. Q which describes the dynamics of such a limit has an infinite slope. In 
this section we will compare numerically trajectories when = and when ^ 0. 
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t(MCS) 

Figure 3. Sample trajectories obtained on a heterogeneous track using the model 
with = for various values of the reduced force f /T. The lower the trajectory the 
higher the opposing force. Here we took with equal probability a = 5,u> = 1,lu' = 2 
and a = 0.2, w = l,u/ = 1. In both cases Ae = and e Afl ^ T = 500. In the inset the 
two largest values of the force are presneted in more detail. 

Using Eq. ^2 the infinite slope of the effective energy landscape implies that 
/i = oo > 2, so that the dynamics is diffusion with drift. The linear drift is illustrated 
clearly in Fig. 01 where trajectories of the strongly biased model are shown for different 
forces. Note that even for very large forces the displacement of the motor as a function 
of time grows linearly. It can be shown that the dwell time distribution on the track 
decays exponentially V(r) ~ e~ T l T * in contrast to the power law distribution expected 
when w£~ ^ 0. To observe any significant pausing clearly one must have / ~ Afi (see 
inset of Fig. EJ). In this limit, however, > w^, and backwards motion cannot be 
neglected (see Eq. IT2^) . 

In Fig. 13] we show the very different behavior of simulations where we take the 
back hopping to be non-zero. Now we do not neglect uo and the backward hopping 
rate is nonzero! For forces near the stall force the displacement of the motor 
as a function of time seems to saturate even for a trajectory generated by a single 
numerical experiment for a particular sequence. Such a behavior is consistent with that 
expected from a sublinear displacement of the motor. Note, however, that with the 
exception of the largest force, all curves on long enough time scales are expected to 
yield, after an average over many thermal realizations, an asymptotically linear curve 
(see the corresponding values of \i presented in the figure). However, even for these 
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Figure 4. Trajectories obtained using the model of Eq. El with = (a + uj)e^ 2T 
and = (ae A,1 / T + uj)e~ Ae / T ~^ 2T . The lower the trajectory the higher the 
resisting force. Here we took with equal probability a = 5, to = I, to' = 2 and 
a = 0.2, lu = l,ui' = 1 along the heterogeneous track. In both cases Ae = and 
gA^/T _ 500. Anomalous displacement occur for f/T > 2.47 while anomalous diffusion 
occurs for f/T > 2.38. The stall force is f s = 2.62. 

curves, pausing is pronounced despite the fact that for small resisting forces the motor 
rarely moves backwards. 

4. Effect of heterogeneous fuel concentrations for a motor moving along a 
heterogeneous substrate 

For some molecular motors the type of fuel which is used to move depends on the specific 
site along the track. For example, in the case of RNA polymerase, which produces 
messenger RNA, the energy from the hydrolysis of the specific NTP which is added to 
the mRNA chain is used for motion. While a random forcing energy landscape would 
exist even if the chemical energy released from every NTP were the same (see Sec. 2), 
the different chemical energies enhance the variance of the slopes, V, of the random 
forcing energy landscape (due to the different chemical potentials of the nucleotides in 
the transcript). This in turn (see Eq. HT|l lowers the value of /i as compared to the case 
of equal chemical energies. 

As suggested in j2J the variance V could be further increased by increasing the 
concentration difference between the different NTPs in the solution. In an extreme 
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situation, where one of the NTPs is completely removed from the solution, the motor 
will stall at specific locations where the NTP is needed. This trick is used in experiments 
to synchronize and control the motion of RNA polymerases [HDj. In this section we 
illustrate using the simple model, Eq. the effects of changing NTP concentrations 
on motor motion. Although we work within a "minimal model" we expect the same 
effects for more complicated models of molecular motors with many internal states (for 
an example of such a model for RNA polymerase see |31j). 

To this end, we studied motors which can use two kinds of fuels depending on its 
location along the track. We hold the concentration of one fuel fixed and lower the 
concentration of the other by reducing the chemical potential associated with it. For 
reference we also study the case when the fuels have equal chemical potentials. Fig. 
El displays results of numerical simulations of the model defined by Eq. |T2|) . Similar 
results were obtained with the more general model of Eq. We show simulations with 
f /T = and when for one fuel e Af */ T = 500 and for the other e A/X//T = 500, 50 or 5. From 
a generalization of Eq. |2l we see that this latter variation corresponds to a change of two 
orders of magnitude in the concentration of the second fuel ^3] . As can be seen from 
Fig. El the effect of reducing one of the chemical potentials is to slow down the motion. 
However, note that for all fuel concentrations the motor displacement as a function of 
time is linear with these parameter values. Since at zero applied force one expects the 
motor always to be biased preferentially in a given direction which is independent of the 
monomer this behavior will hold even for larger differences in concentration. 

The change in concentration of one of the fuels can make the regime of anomalous 
dynamics much larger. In Fig. |U]we show the result of applying a force which opposes 
the motion of the motors on the dynamics. As can be easily seen from the figure, 
the velocity of the motors is, of course, lower in comparison to that when no force 
is applied. However, note that the motion of the curve with the largest difference in 
chemical potential (barely visible at the bottom) is very different. The motor almost 
immediately stalls after it starts moving. This example illustrates clearly how changing 
the concentrations of the different fuels can make the regime of anomalous dynamics easy 
to access at lower forces. In fact, averaged over many thermal and sequence realizations, 
the curve with the largest difference in chemical potentials shows a displacement of the 
motor which grows sublinearly in time, indicating that /i < 1. 

5. Heterogeneous fuels on periodic tracks 

In this section we consider a molecular motor which is moving along a periodic track 
in a solution which contains more than one type of molecule which can supply it with 
chemical energy. The situation arises for kinesin in the presence of both ATP and, say, 
GTP. It is known that, while less efficient, alternative NTP molecules can also be used by 
kinesin to move along the track [231 EH EH] • We consider, for simplicity, a solution with 
two types of chemical fuels within the simple two state model (The analysis presented 
can easily be generalized to include additional fuel types). In addition, we generalize 
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Figure 5. Motor trajectories obtained using the model denned by Eq. 1|12|) . Here, 
substrate heterogeneity is incorporated by taking with equal probability {a = 5,c<j = 
1, uj' = 2} and {a = 0.2, u> = 1, u' = 1}. The corresponding values of e A/ ^ T for the two 
different fuels are denoted in the box. In both cases Ae/T = and f/T = 0. Lower 
trajectories correspond to a larger difference in the chemical potentials associated with 
the fuels. 

our model to treat two separate cases. In the first we assume that the internal states 
of the motor (i.e., the a and b sites in Fig. UJ) are independent of the fuel used so that 
the chemical fuels are only used to move between states with fuel- dependent potential 
differences driving the changes. This situation arises when the fuel from the motor is 
used and released so quickly that the motor is unbound to the fuel in the "excited" 
internal state. In this case the fuel binding to the motor does not define an internal 
state. In the second case we study we allow for additional internal states of the motor, 
depending on which type of chemical fuel is bound to it. This situation arises when 
internal "excited" states include a fuel bound to the motor. Since the motor bound to the 
two fuels defines two distinct internal states, a direct thermal transition between them 
is therefore not possible. While both cases discussed below involve parallel pathways 
for transitions across a monomer, they are distinct in the type of internal states. 

5.1. Case I: Internal states of the motors are independent of the chemical fuel 

This case amounts to a straightforward generalization of the rates in Eq. Q to allow 
for multiple fuels. As mentioned above, we assume the internal states of the motor are 
independent of the chemical energy and that chemical energy is only used to assist the 
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Figure 6. Motor trajectories obtained using the model Eq. I|12f> with an opposing 
force. We again took with equal probability {a = 5, u> = 1, u>' = 2} and {a = 0.2, uj = 
1,uj' = 1}, with the corresponding values of e Afl ^ T for the two different fuels as 
indicated in the figure. In both cases Ae/T = and f/T — 1. Lower trajectories 
correspond to a larger difference in the chemical potential differences. 

transitions. With two different chemical fuels, the generalized rates entering Eqs. Eland 
Eq. H]are now 

= ( ai e A ^ T + a 2 e A ^ T + uj)e~ A ^ T ^l 2T 
= («i + a 2 + uj)e f/2T 
w ^ = ( a [e A ^ T + a' 2 e A ^ T + a/)e" A£/T+//2T (13) 
= (a[ + a' 2 + uj')e~ f/2T . 

There are now three parallel paths, two assisted by chemical energy, one thermal. The 
subscript 1 or 2 refers to whether fuel "1" or "2" is being used for the transition. 
Standard relations for the chemical potential difference ^3] imply that e AMi//T grows 
linearly with the concentration of fuel "i" (see Eq. EJ). Because the extra channel is 
present, when one of the chemical potential differences is set to zero the model reduces 
to the single motor model but with modified rates as compared to a motor in a solution 
where the fuel is completely absent. For example, if Afi 2 = the model can be written 
in terms of Eq. ^with the identifications u — > uj + a 2 and uj' — > uj' + a' 2 . Note also 
that, in contrast to disordered tracks where the motion of the motor is described by 
a random walker moving on a random force landscape, here the energy landscape is 
periodic except for a well defined tilt given by Eq. 
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With the help of Eq. |SJ it is straightforward to verify that for this model the 
velocity takes the form 



where [x±] and [xz] are the concentrations of fuel 1 and 2 respectively and 
A, B, B', C, D, D' are complicated functions of the the coefficients in Eq. ED The 
velocity thus varies as ratio of two polynomials, each depending linearly on the 
concentration of both fuels. Consider, for example, the velocity in an experiment where 
one of the fuels is held at constant concentration and the concentration of the other 
is varied. The external force is set to / = 0; (varying / did not affect the qualitative 
features discussed below). As shown in Fig. the average velocity smoothly crosses 
over from a small value to a larger value in a sigmoidal fashion as the concentration of 
more energy rich fuel ("fuel 1") is increased. Note that A/i 1( /T varies logarithmically 
with the concentration of fuel number 1 [T^j. Thus, the concentration of fuel "1" varies 
over many orders of magnitude for the range of A/xi/T shown in Fig. In typical 
experiments, only a small portion of this crossover may be visible. For small A/xi/T, 
motor movement is controlled by fuel "2" while for large A/xi/T it is controlled by fuel 
"1". An analogous plot for motor velocity vs. Ajiii/T when fuel 2 is absent entirely is 
shown for reference in Fig. |H1 As expected, the velocity no longer exhibits a sigmoidal 
crossover between two regimes. 

5.2. Case II: Internal states of the motors are coupled to the chemical fuel 

Next, we consider a different model incorporating two fuels. We now assume that the 
type of fuel molecule used to move the motor determines the entire chemical cycle which 
leads the motor to move across one monomer. An example might be the n = 2 motor 
landscape shown in Fig. Qwhere one step (e.g.; a site — > b site) involves the fuel molecule 
binding to the motor. In this case the entire sequence of transitions would be dictated 
by the fuel that is used. Thus, the choice of rates (either in the forward or backward 
direction) would be dictated by the initial step which is chosen at random, depending 
on the type of fuel utilized. 

Since the state of the motor is directly coupled to the chemical fuel, a pure thermal 
transition between the states is not possible (unless it involves moving across the 
monomer through parallel pathways not related to the motor, an effect which is ignored 
here). The new feature is that the motor can move across a monomer by choosing one 
of two distinct chemical pathways. This is in contrast to the model defined by Eq. fl~3l 
where two distinct chemical pathways exist for passing between the internal states of the 
motor, but where the transition across the monomer can occur via a mixture different 
chemical (or thermal) pathways. 

An example for such a model, where each fuel is modeled by a distinct channel 
which is a special case of the two state model defined by ((TJ) is given by: 



v = 



A + B[ Xl } + B'[x 2 ] 
C+D[xi] +D'[x 2 ) ' 



(14) 



= ae 



A M i/T-//2T 



a 
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Figure 7. The velocity as a function of the chemical potential difference A^i/T 
of fuel "1". The velocity is plotted in units where is size of a monomer (i.e., the 
lattice constant 2ao in Fig. ^) is set to be one and arbitrary time units. We used 
a "nearly one way model" based on Eq. ^] where the rates were chosen to be 
a\ = 5,Q'2 = 2, a[ = 0.01, a' 2 = 0.02, uj = 0.01, uj' = 0.2 in arbitrary units of inverse 
time and A^/T = 5, / = 0, Ae = 0. For reference the velocity of with only fuel of 
type "2" with the same value of A/X2/T is v = 0.195. 



= ae^ 2T 

w^ = cue f/2T (15) 



for the first channel and 

u7 = le A, 2 /T-f/2T 

= -ye^ 2T 

u ^ = i/e //2T (16) 
u - = v e-f' 2T . 

for the second channel. The fuel concentrations enter through the chemical potential 
differences A/ix and A/i 2 (see Eq. |2J). The relative fuel abundances therefore control the 
ratio of the rates and . The model will be realized physically in a motor which 
binds a fuel and uses its chemical energy in the transitions or . The transitions 
and involve the release of the corresponding fuel. Here for simplicity we have 
set the energy difference between the states to be zero and neglected parallel thermal 
channels. We do not expect the qualitative results described below to be affected by 
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Figure 8. The velocity in arbitrary time units as a function of the chemical potential 
difference A/ii/T where the other chemical channel associated with fuel "2" is closed. 
The velocity is plotted in units where is size of a monomer is set to be one. The rates 
were chosen to be ot\ — 5, a 2 = 0,a[ = 0.01, a' 2 = 0.0, a; = 0.01, w' = 0.2,/ = in 
arbitrary units of inverse time and Ae = 0. 



such complications. 

The velocity of the model can be calculated in a straightforward manner and is 
found to be 

v = (v + ^rx^rv - KjjC) + or + ^rx vv - <<) 

(V + U V)( W V + W D + ( U V + K)( w a + w l) + (?C + u l)( w V + W D 

Note that as in Case I, the velocity behaves as a ratio of two polynomials which are 
linear in the e A/ ^/ T , and hence with each of the fuel concentrations (similar to Eq. IT^j) . 
Also similar is the fact that the presence of a second channel alters the velocity even 
when one of the chemical potential difference is zero. Again, due to the presence of the 
second channel the velocity is lower than that of the motor in the presence of a single 
fuel. We do not present plots of the resulting velocity as the concentration of one of 
the fuels is varied since the qualitative features are similar to those presented in the 
previous subsection. 

Finally, we comment that more complicated scenarios (for example, the existence of 
additional parallel pathways, or more general n-state models) might change the explicit 
dependence on the concentrations of the fuels. However, we expect a general form of a 
ratio of two polynomials in the concentrations of the fuels even in much more complicated 
scenarios. Indeed, with a specific experiment in mind and some structural information 
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on the motor one might be able to use such experiments, coupled with an analysis 
similar to that presented above, to deduce the number of steps in the chemical cycle 
which depend explicitly on the fuel concentration. For example, we have considered 
models where the internal states are independent of the chemical fuel and found that 
the general velocity can be a ratio of polynomials of a degree which is related to the 
number of steps which depend on the chemical fuel. However, the detailed results were 
very dependent on the choice of rates made. 
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